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I. INTRODUCTION 

Onsager’s celebrated paper [l| on the effect of shape on the interaction between hard particles has perhaps been 
the most influential contribution to colloidal sciences of the last century There, entropic forces alone were first 
recognized as capable of inducing a structural ordering transition with no involvement of whatever cohesion force may 
be present. The typical prototype of such an ordering transition remains indeed the isotropic-to-nematic transition 
predicted in Q for an assembly of slender hard rods as their number density is increased beyond a critical value 
(falling within a narrow gap of phase coexistence). As paradoxical as it may appear at a superficial glance, such an 
ordering transition is duly accompanied by an increase in entropy, since the loss in orientational disorder attached to 
the rods’ alignment is outbalanced by the gain in translational disorder made possible by the increase in the volume 
available for the particles’ centers of mass [I, Q- The conjugated counterpart of this volume is the excluded volume. 

The excluded volume of two rigid bodies is the volume in space that any one point in one body cannot access by 
the very presence of the other body. This definition is delusively simple as it conceals a formidable mathematical 
task which can seldom be accomplished in an exact analytic form 0 Of course, there are exceptions to this general 
statement but they are very few. Noticeable among these are the excluded volume of circular cylinders [l|, sphero- 
cylinder [ 6 [, sphero-platelets Q, and sphero-zonotopes 00 

Despite its technical difficulties, the excluded volume remains a key ingredient of both Onsager’s original theory and 
its most recent extensions. In all of these, the per-particle free energy F of an assembly of hard bodies (appropriately 
made dimensionless) is a functional of the single-body local density g. A number of papers have interpreted Onsager’s 
original theory in the light of the modern density functional theories; here we refer the reader to the most recent 
review on the subject jl^, which is mostly concerned with hard-body systems that exhibit liquid crystalline phaseslH 
F[g] differs from the free-energy functional for an ideal gas by the addition of an excess free energy Fex[^], which 
characterizes the interactions of anisometric particles. In general, Tex[^] is not known explicitly, but it can always 
be expressed as a power series in the total number density po, which is often called the virial expansion. The first 
non-trivial term of such an expansion is poB2[p], where the functional B2 is the second virial coefficient, which is 
nothing but the ensemble average of the excluded volume, 

^ 2 [q]'= 7 ; [ Ve{uj,uj')g{uj)g{uj')dujduj'. ( 1 ) 

In ([ 1 ]), U is the orientational manifold^ which describes all possible orientations of a particle in the system and 14 (cj, uj') 
is the excluded volume of two particles with orientations uj and cj', respectively. Higher powers of po bear higher virial 
coefficients B^, which however are even more difficult to compute than B2. 

Onsager [l| remarkably estimated that for rods sufficiently slender B2 actually prevails over all other B^’s. This 
makes Onsager’s theory virtually exact, as was also subsequently confirmed directly by numerical computations 
mui. Nevertheless, even when the second virial coefficient B2[g] cannot be proved to be dominant, it remains a 
viable approximation to Fex[p] in establishing, at least qualitatively, the variety of possible equilibrium phases in a 
hard-body system and the entropy-driven transitions between them. To this end, explicit formulas for the excluded 
volume of rigid bodies are to be especially treasured. 

This is the motivation for our study. Our objective is to express ^2}^ the excluded volume for two rigid 

bodies, and ^2, in terms of shape functionals depending solely on the individual bodies and ^2- We shall 
accomplish this task for bodies both convex and cylindrically symmetric, for which Ve{^i,^2} can be given with no 
loss in generality as the sum of a series of Legendre polynomials Pn , 

c» 

Ve{^l, ^2} = 'F BnPn{mi ■ m2), ( 2 ) 

n=0 

where mi and m2 are unit vectors along the symmetry axes of and ^2, respectively!! The shape functionals 
involved in our explicit representation will be natural extensions of the classical functionals on which was largely 
based the celebrated Brunn-Minkowski theory of convex bodies]! The major advantage of the method proposed here 
is the explicit computability of such extended Minkowski functionals, which makes our representation formula directly 
applicable to bodies and ^2 not necessarily congruent, possibly representing particles of different species. 


^ We learn from Q that Viellard-Baron, who took an early interest in this problem “was reportedly greatly disturbed by the difficulties 
he encountered.” 

^ Isihara @ is often credited with having provided an explicit formula for the excluded volume of ellipsoids of revolution. In Sec. EH 
below, we shall discuss this case in some detail. 

^ A general reference for simple liquids is still the classical book now enriched by an addition on complex fluids. 

^ Following Isihara Q, we denote by Bn the Legendre coefficients of W, though often in more recent literature this symbol is used to 
designate the virial coeffients, here denoted as Bn- 

^ Besides the original sources [ij, [3 5 the general books EE US are highly recommended. We also collected a number of relevant results 
phrased in the same mathematical language employed here in Appendix A to our earlier study on this subject El- Finally, a different 
but equivalent approach is presented in [l^ . 
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The paper is organized as follows. In Sec. [Ill we set the scene for our development by showing that the Legendre 
coefficients of the representation formula m can be expressed as appropriate anisotropic volume averages. Sec¬ 
tion M is devoted to the coefficient Bi of the first Legendre polynomial Pi (mi • m2) = mi • m2 in ([ 2 ]). We attach 
a special meaning to this, as it represents the dipolar contribution to ^2} which would possibly arise from 

tapered, cylindrically symmetric, convex bodies, if only one could unambiguously assign a shape dipole to them. The 
somewhat surprising conclusion will be that Pi vanishes identically on this class of bodies, making the very notion of 
shape dipole void, despite its intuitive appeal. Section llVl is concerned with the extended Minkowski functionals, in 
terms of which, once evaluated on the bodies ^1 and ^2, we can write in closed form all coefficients B^ in ([ 2 ]). An 
explicit application of our method is illustrated in Sec. El where we evaluate the extended Minkowski functionals for 
a generic circular cone and validate our evaluations through a direct computation of the coefficients B^ made pos¬ 
sible by an independent shape-reconstruction algorithm, appropriately modified to tackle efficiently the cone’s sharp 
ridge. Likewise, in Sec. [Vll we determine the extended Minkowski functionals for a spheroid, that is, an ellipsoid of 
revolution, either prolate or oblate. In Sec. Em we collect the main conclusions of our work, looking back afresh to 
some of them, also in the light of possible future developments that they may suggest. 

We shall endeavor to make our presentation as free as possible from unwanted technical details that might obscure 
both the outcomes of our study and the strategy adopted to obtain them. To provide, however, the interested reader 
with enough information to appreciate the mathematical infrastructure underlining this paper, we collect in two 
closing appendices the details of both the mathematical theory and the shape-reconstruction algorithm. 


II. ANISOTROPIC VOLUME AVERAGES 

It was proved by Mulder [ 8 [ that the excluded volume of ^2} of two bodies, and ^2, be they convex or 

not, can be expressed as 


K{^1,^2} = I"[^1+^2*], (3) 

where V is the volume functional, ^2 the central inverse (relative to a specified origin o) of the body ^2, and 
+ denotes the Minkowski addition (to the definition of which concurs the origin o)E Letting both and ^2 bo 
cylindrically symmetric bodies with axes mi and m2, respectively, since ^2} is an isotropic scalar-valued 

function, by a theorem of CauchyQ we can say that ^2} is a function (still denoted as) W of the inner product 

mi • m2. Setting mi • m2 = cos'd, the function K(cos'd) can be expanded as the sum of a series of Legendre 
polynomials (see, for example. Secs. 18.2 and 18.3 of | 2 l|): 


00 

14 (cOS??) = V BnPnjcOS'd), 

n=0 


(4) 


where 

277 + 1 

Bn'.= - / K(cos'd)Pn(cos'd) sin'dd'd ( 5 ) 

2 Jo 

are the Legendre eoeffieients of Vq. We record for future use a few basic properties of the orthogonal polynomials Pn 
(see, in particular. Secs. 18 . 6.1 of M and 8 . 917.1 of My- 

Pn{-X) = i-irPnix), P„(l) = l, |P„ + +1. (6) 

There is another way of expressing the coefficients Bn^ which we find illuminating. Consider the average 

(PnVe) [^1,^2] := {Pn{mi • m2)Ve{mi •m 2))^2 (^) 

computed for fixed ^1 over all possible replicas of ^2 obtained by rotating arbitrarily ^2 in space. By the cylindrical 
symmetry of ^2, the average 0 also acquires the equivalent form 

(PnVe) = {Pn{mi ■ m2)Ve{mi ■ m 2 )), ( 8 ) 


® We shall often call 0 Mulder’s identity. The reader is referred to the primer on the Brunn-Minkowski theory of convex bodies in 
Appendix A of [l^. A short recapitulation of this theory is also given in Appendix ixn below to make our paper self-contained. 

^ See, for example, Sec. 113.1 of [20l| . 
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where, for any function /(e) defined on the unit sphere 

(/)e := / / fie)daie) (9) 

47r yg2 

and da{e) denotes the area element with unit normal e. Representing m 2 in polar spherical coordinates with polar 
axis mi and combining ^ and ([5]), we readily arrive at 

(PnVe) [^ 1 ,^ 2 ] = ^ / Pn (cos'd) Re (cos'd) sin 'dd'd = ^ Pn> (10) 

2 do 2n + 1 

Since both functions Vq and Pn are symmetric under the exchange of mi and m 2 , the average (PnK) [^ 1 , =^ 2 ] is also 
symmetric under the exchange of bodies ^1 and ^ 2 ' 

(P^Re) [^ 1 ,^ 2 ] = {PnVe) [^2,^l]. (11) 

Equation © allows us to express the Legendre coefficients Bn of the excluded volume of two cylindrically symmetric 
bodies in a way directly related to the anisotropic averages of the volume of a Minkowski sum. Combining (|TU1). ©, 
and ©, we readily see that 

Bn = (2n + 1) (P„(mi ■ m 2 )V[^i + =^ 2])<^2 = ■ m* 2 )V[^i + 

= (2n + 1)(-1)” (P„(mi ■ m* 2 )V[£gi + 

where m^ = —m 2 is the symmetry axis of the central inverse ^2 ^^0 use has been made of (|6|) and the fact 

that averaging over ^2 is just the same as averaging over Thus, to obtain all coefficients in (|4]), we need to 
learn how to compute the anisotropic volume averages 

(PnV) [^ 1 ,^ 2 ] := (PnR[^l +^ 2 ])^, , (13) 

as then ([T2]) would simply reduce to 

Bn = (2n + !)(-!)" {PnV) [^ 1 , ^ 2 *], (14) 

which obeys the same symmetry relation as in (pT]) . Equation (p!4|) is the basic building block of our development. 

Although da is as general as © for cylindrically symmetric bodies, this paper will solely be concerned with the 
excluded volume of convex cylindrically symmetric bodies. Eor n = 0, the average in m becomes isotropic as Pq = 1 
and its expression has long been know for generic convex bodies 0 

(y) [^ 1 ,^ 2 ] = V[^i] + V[^2] + T {M[^^]S[3§2] + M[^2\SWi\) , (15) 

47r 

where M is the total mean curvature functional in (|A12aP and S is the surface area functional in (|A12bp . Since both 
M[^] and S[^] are invaraint under central inversion of it follows from (p!T|) and ([T5|) that 

^0 = (4^) [^ 1 ,^ 2 ]. (16) 

Here our challenge is to extend the neat classical formula (p!5]) for the isotropic average of the volume of the Minkowski 
sum of convex bodies to the anisotropic averages needed in (p!T|) . This will be achieved in the two following sections 
with the aid of appropriate extensions of the classical Minkowski functionals M and S. We anticipate that they are 
invariant under central body inversion like the classical Minkowski functionals, so that, in complete analogy with m 
and (HD, we shall be able to express the excluded volume R{^i,^ 2 } of cylindrically symmetric bodies and ^2 
in terms of functionals evaluated separately on and ^ 2 - 

As recalled in Appendix El there is no loss in generality in limiting attention to the class of convex bodies 
with smooth boundaries and strictly positive principal curvatures, as is dense in the whole class of convex 
bodies (see Appendix lA ip . Thus, our strategy will be to compute first the anisotropic volume averages in and 
then extend them by continuity to the whole of JV. In the following section, we shall first accomplish our task for 
(PiR) [=^ 1 ,^ 2 ]; this will lead us to conclude that Pi = 0, a general result of some import. In Sec. [IVl we shall 
compute {PnV) [^ 1 ,^ 2 ] for all n ^ 2 and arrive at the expected general explicit formula for all P^’s. 


® A derivation of (|T5) i can be found in 0. Moreover, Kihara [^. credits Isihara and Isihara and Hayashida |26|. |27 
proved CD, although he also seems aware that a proof had already been contained in the classical work of Minkowski \l5 


for having 
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III. NO SHAPE DIPOLES 

Here our task is to compute Bi. To this end we remark that 

{Pn{mi • m2))^^ = {Pn{mi • m2))^^ = 0, (17a) 

{Pn{mi • m2)V[^2])^^ = V[^2] {Pn{mi • m 2 ))6^^ =0, n ^ 1, (17b) 

the former following from (j9]) and the orthogonality of Legendre polynomials, and the latter also from the invariance 
of the volume functional under rotations. Then we represent m 2 in a Cartesian frame (e^,, ey^e^) fixed in ^ 1 . Letting 
Bz = ly, where ly is the outer unit normal to at a selected point on and choosing By orthogonal to the plane 

(mi, I/), we have that 

mi = sin'diCa, + cos'll L/, ( 18 a) 

m2 = cos (p sin '^263, + sin 0 sin '^26^ + cos '^2^', ( 18 b) 

the latter of which represents all possible orientations of m 2 , for given 'di and '^ 2 , the angles that mi and m 2 make 
with ly (see Fig. [T]). An easy, but important consequence of ([18]) is that 


mi 



m 2 


FIG. 1. Sketch representing the unit vectors i/, mi, and m 2 . With ly and mi fixed, m 2 as represented by ()18bl) describes a 
cone around u in the first step of the averaging process described in the text. 


mi • m 2 = sin 'di sin 1^2 cos (j) + cos 'di cos 1^2 

= sin 'di sin 1^2 cos <p + (mi • ly) {m 2 • jy ). 

Now, using also (13, we can derive from (jA14j) the following expressiorU 

{PiV) [S§i,S§2] = ^ ^ ■ ri){u ■ rni)da{u) + {{u ■ r2){u ■ m2))^ ^\i)^ da{u^ 

+ i {{m 2 -v) {pi^ + p^ 2 ^yj ^ ■ ri){mi -u) {p^^^ + p 2 ^) da{u) ( 20 ) 

+ ^ ((i^ • r2){m2 • v) + P2^))^ 

which results from computing the average over ^2 in two separate steps: first averaging over the angle (j) in ([T9|) 
which ranges in [0, 27r] and then averaging formally over i/, meant as the outward unit normal to =^ 2 , which ranges 
over §^. If the former average is taken over the process in which, with u and mi fixed, m 2 is seen to describe a cone 
around ly (see Fig. [T]), the latter is nothing but the average over the independent process in which all different points 

of d ^2 come to be associated with one and the same fixed normal ly. As in (|A14p . also in (Eni and denote the 

(2) (2) 

principal radii of curvature of d^i and pi and p 2 denote the principal radii of curvature of d^ 2 ] correspondingly, 
1 ^( 1 ) = (p^pp^P)-'^ and ^( 2 ) = (pPpP)- are the Gaussian curvatures of d^i and d ^2 and ri and r 2 are the radial 
mappings of ^1 and ^2 (see Appendix I A II for more details). 

Now, with the aid of the theory recalled in Appendix I A 11 we compute the new shape functionals featuring in (|2Q]) . 
It readily follows from (|A8p that for any body ^ G 

/ ^ da{iy) — / m’nda{n)= / divmdu = 0, (21) 

K Jq^ 


^ Unlike Mulder’s identity which is valid for general bodies, equation lIXull . which is indeed one basic ingredient of our theory, has 
only been established for convex bodies. 
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where use has also been made of the classical divergence theorem (and the fact that m can be extended to the whole 
space as a uniform field). Likewise, (|A 9 p and (jASp imply that 


[ {m ■ v'){pi ^ p2)da{u) = [ {m ■ divsn da{u) 

y§2 K 

= / (m • n) divs n (ia(n) = / divg m (ia(n) = 0, 
Jd^ Jd^ 


( 22 ) 


where use has also been made of the surface divergence theorem recalled in (jAlip . Combining (|2T]) and (|22]) . we 
obtain from (j 20 )) that (PiV) [^1,^2] vanishes identically for all and ^2, and so, by ([H]), 


5i = -3(Piy)[^i,^2*]=0. (23) 

Equation (| 23 |) says that for cylindrically symmetric bodies, and ^2, the excluded volume in (j 4 ]) does not 
contain any dipolar contribution, no matter how tethered and ^2 can be, suggesting that no shape dipole can 
be associated with them. It was already argued in that a shape dipole cannot be unambiguously assigned to a 
body Equation (| 23 |) shows that no matter how we endeavor to assign a shape dipole to ^ it plays no role in 
the hard-particle interactions governed by the excluded volume. Of course, polarity effects are also expected to be 
seen in these interactions. Eor example, it was proved in [ 28 | that the excluded volume of two congruent cylindrically 
symmetric convex bodies is minimized when the bodies are in the antiparallel configuration, where m2 = —mi. 
Such polar effects, however, cannot involve shape dipoles: as shown in [^, they start being manifested through the 
shape octupole that features in (| 1 ]) through the coefficient P3. This and all higher order Legendre coefficients will be 
computed in the following section. 


IV. EXTENDED MINKOWSKI FUNCTIONALS 


Computing the anisotropic volume averages (PnV) [^1, ^2] for n ^ 2 is technically more complicated than comput¬ 
ing {PiV) [^1, =^2], although conceptually this task is not much different from that just accomplished in the preceding 
section. As shown in Appendix lA 21 this computation led quite naturally to the introduction of a number of shape 
functionals that extend the classical Minkowski functionals M and S. They are defined for all n ^ 2 as follows: 


Mr,m := 

/ Pn{m ■ n)Hda{n), 

Jd^ 

(24a) 

:= 

/ {n ' x)Pn{m ’ n)Kda{n)^ 

Jd^ 

(24b) 

:= 

[ [1 - (m ■ n)2]i((Ti - ■ n)da{n), 

Jd^ 

(24c) 

Sum ■= 

/ Pn{m ' n)da{n), 

Jd3S 

(24d) 

S'nim ■■ = 

/ {n ' x)Pn{m ' n)Hda{n), 

Jd3S 

(24e) 

s':m ■■= 

[ (n • x)[l - (m • n)^]^(cri - (T 2 )P^^_ll\m • n)da{n). 

Jd3S 

(24f) 


We shall often refer to them as the extended Minkowski functionalsThey give (PnV) [^1, ^2] the following concise, 
explicit representation: 


(PnV) {M^[m]Sn[^2] + + T (M„[i^i]5;[^2] + M„[^2]5;[i^i]) 

i27r OTT 


1 (n —2)!(n + 2)! 


Gtt 


( 4 n!) 


!i2 


(M"[^i]5"[^2] + M"[^2]5"[=^i]) . 


(25) 


Strictly speaking, in Appendix lA 21 we arrived at (|M|) through the representation via radial mapping of the convex 
bodies in the special class . However, the extended Minkowski functionals can also be extended by continuity to 


10 


More shortly, also as the extended M and S functionals. 
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the whole of . Moreover, as clearly shown by (|2^ . their definition actually applies to any cylindrically symmetric 
body, be it convex or not. The extended M and S functionals are invariant under rotations. Their behavior under 
translations is further discussed in Appendix IA51 

Since the extended Minkowski functionals for a body ^ are invariant under central inversion of ^ (see Ap¬ 
pendix [AT]), it follows from (|25]) that (PnV) = (PnV) [^ 1 ,^ 2 ], and so equation (fT4|) becomes 

Br, = (2n + (PnV) [^ 1 , ^ 2 ], (26) 

which by ([25]) expresses the Legendre coefficients of 14 in (|4|) in terms of shape functionals evaluated on the individual 
bodies and ^ 2 - Formula (j26|) will be applied in the two following sections to special classes of bodies, namely, 
circular cones and ellipsoids of revolution. 

As shown in Appendix I A 61 the functionals M^', M^, and are not independent of one another. If is certainly 
related to through 

= vn^2, ( 27 ) 

77/ —n ^ 

for all cylindrically symmetric convex bodies we expect the relation 

= - (n-l)(„ + 2) “”W 

to be valid at least for both classes of bodies studied in detail in this paper, having checked it by direct inspection for 
a large number of indiceslHI Whenever (|28|) applies, the anisotropic volume averages {PnV) [^ 1 ,^ 2 ] in (|2S|) take on 
a much simpler form, 

{PnV) [m,^2] = T {Mn[m]An[^2] + Mn[^2]An[m]) , (29) 

OTT 

where 

1 _l_ 1 

:= V n ^ 2 . (30) 

In particular, for two congruent bodies, ~ ^2 ~ =^0 by ([26|), Bn can be given the following factorized expression, 

^ (2!^±^KzDlMn[^]Anm, (31) 

OTT 

which we shall assume to be valid in the following (and will be very convenient in our development below) 0 


V. CIRCULAR CONES 


We denote by a circular cone with semi-amplitude a G [0, ^], radius R, and height h, both related through 
(|A32p to the slant height L (see Fig.[2|). It is a simple matter to show that the classical Minkowski functionals for 



FIG. 2. (Color online) A circular cone with vertex in the origin o, semi-amplitude o, radius R, height /i, and slant height L. 


Of course, we are aware that this can by no means be considered as a proof of CHli, which remains for us a conjecture, though with a 
high likelihood of being true. 

Meaning that and ^2 are images of ^ under the action of the full orthogonal group 0(3). 

In the language of [2^ and [^, once combined with di, EDI would be called a convolution decomposition (or simply a deconvolution) 
of the excluded volume. 
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take the explicit forms (see also (A61) and (A62) of (ill), 




ttL 


cos Of + 




sin a 


= ttL^ sin(a(l + sin a), 


V[^^] = -ttL^ cosasin^ a. 

o 


(32a) 

(32b) 

(32c) 


As follows easily from (jA34p . the Gaussian curvature K vanishes identically on all smooth components of 
Moreover, the contribution of the vertex o to all the integrals in (|M|) vanishes, as can be seen by replacing o with a 
fitting spherical cap of radius e (whose area surface scales like e^) and then taking the limit as 5 ^ 0^, in complete 
analogy to the method used in Appendix IA4al to compute the extended Minkowski functionals on a circular ridge 
93. The formulae (jA4Q[) obtained there for a 93 can be directly applied here to the rim of the cone’s base by simply 
setting 0i = ^ — a and 62 = 71 . Use of (|A32p finally leads us to 








ttL yP, 
-27rL 


(sin (a) cos (a + sin (a / Pn{cos'd)d'd , 

h-^ J 

r 

/ cos('d + OL)Pn (cos'd) sin 'dd'd, 
df-a 

—ttL (sin (a) cos^ (a — sin (a J Pn-^2 '^) 'dd'd^ 

71 sin Of [Pn(sin(a) + (—1)’^ sin a] , 

—7rl/^sin(a / cos('d + (a)Pn(cos'd)(i'd, 

J^-a 

pTT 

—ttL^ sin a / cos('d + (cos 16) sin^ 'dd'd, 

J ^ —a 


(33a) 

(33b) 

(33c) 

(33d) 

(33e) 

(33f) 


for all n ^ 2. Inserting (IMl) in (IMll . we obtain explicit, analytic formulae for the Legendre coefficients Bn of the 
excluded volume of two congruent circular cones, and ‘^2 ^ which for completeness are recorded in (IA51 1) for the 
first seven indices n ^ 1. They are plotted in Fig. [3] as functions of a. Inserting (11^ in (USD, we also obtain the 
isotropic average Pq in ([T 6 |) . which is plotted in Fig. |4] with two possible normalizations, relative to the volume Vc 
of each cone delivered by (|32cp in Fig. HJ^a), and relative to in Fig. (U^b). The even-indexed coefficients are 
mostly negative, indicating by m a tendency for the corresponding terms in the sum to minimize Vq for either 
'd = 0 or'd = TT, irrespectively. On the contrary, the odd-indexed coefficients are mostly positive (apart from P 3 which 
is never negative), indicating a tendency for the corresponding terms in (| 1 ]) to minimize Vq for'd = tt, that is, when 
the cones and ^2 antiparallel configuration, with m 2 = —mi. This suggests that the excluded volume 

of two congruent circular cones is minimized in the antiparallel configuration, as shown by direct computation in [T^ 
in accord with the general minimum property established more recently in . 

The crosses superimposed to the graphs in Fig. [3] represent the values of Bn extracted numerically from the 
volume of the excluded body ^^ 2 }^ region in space that cone ^2 cannot access by the presence of 

cone Determining ^^ 2 ^ is indeed necessary for a direct determination of as the general 

proper geometric definition of the excluded volume of bodies and ^2 is precisely the volume of the excluded 

body •= '^2}] (scc also Hi)- Here ^^ 2 ^ was obtained from the shape- 

reconstruction algorithm outlined in Appendix [Bl Our strategy was completely different from that adopted so far in 
this paper. For a given a, we reconstructed .,^ 2 ^ for a number of values of the angle'd made by the cones’ 

axes mi and m 2 ; we computed numerically the excluded volume I 4 as a function of'd by applying (|A 12 cP to a 
triangulation of '>^ 2 ^ extracted from this function the coefficients Bn through (|5]). To what extent 

the two methods agree, thus granting support to each other, is left to the reader to judge from Fig. [3l Quantitative 
details about both the shape-reconstruction algorithm employed here (including its adaptation to the specific case 
of cones, which with their sharp edge and pointed vertex required special attention) and the way the coefficients Bn 
were computed can be found in Appendix [Bl below. 

Figure [5] shows three graphs representing the excluded volume Ug of and scaled to their common volume Vc 
(given by (j32cp j as a function of the angle'd between their axes. The semi-amplitude a of both cones is taken here to 
be 0^0 = 0.14 TT, for which, as shown in Fig. HI the isotropic average (I 4 ) scaled to Vc takes on its minimum value. The 
graphs in Fig. [5] correspond to the function in ([12]) truncated at n = 3 and n = 9; they are both contrasted against 
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FIG. 3. (Color online) (a) For two congruent circular cones, "ioi and "^ 2 “, with slant height L and semi-amplitude a, the graphs 
of Bn scaled to are plotted against 0^(a^^forn = 2 (solid line), n = 4 (dashed line), and n = 6 (dotted line), according 
to (IA51|) . (b) For the same cones, "ioi and "^2 ^ fhe graphs of Bn scaled to are plotted against 0^(a^^forn=l (thin solid 
line), n = 3 (solid line), n = 5 (dashed line), and n = 7 (dotted line). In both panels, crosses represent the values computed 
numerically on the shape of the excluded body ^^ 2 } reconstructed with the algorithm recalled in Appendix [B] 




FIG. 4. (Color online) (a) The isotropic average Bo as in (flGl) and (ITSl) normalized to the cone’s volume Vc delivered by 
(|32cl) : it attains its minimum at a = 0.147r. (b) Bq normalized to like all other coefficients Bn’s shown in Fig. [Sj it 

attains its maximum at a = 0.477r. In both panels, crosses represent the volumes computed numerically to benchmark the 
shape-reconstruction algorithm described in Appendix m 


the octupolar approximation^ which in [T8| was shown to be rather accurate. While, by construction, the latter takes 
on the exact values of W at both 1 ^ = 0 (parallel cones) and 1 ^ = tt (antiparallel cones), which are 14W and SW, 
respectively, both truncated expansions do not. Actually, as expectedFl the convergence of the series in (|lj) at these 
points is rather slow: for example, a computation with 61 terms was required to obtain 

h = 14,01 and F = 8.153, (34) 

Vc Vq 


Since the expansion in 6) is an approximation in the L^-norm, and not pointwise. 
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FIG. 5. Excluded volume Ve of two congruent circular cones, and ^ 2 , with slant height L and semi-amplitude ao = 0.14 tt 
corresponding to the minimum value of the scaled average (Ve) /Ec, where Vc is the volume of each cone. Two graphs, plotted 
against the angle 0 ^ ^ ^ tt made by the cones’ axes, are delivered by @ truncated at n = 3 (solid line) and n = 9 (dashed 
line). The third graph (dotted line) represents the octupolar approximation proposed in [3, which interpolates the excluded 
volumes of parallel (^ = 0) and antiparallel (^ = tt) conhgurations. 


at '^9 = 0 and i9 = tt, respectively. Thus, if for cones the explicit octupolar approximation of the excluded volume 
could still be a good choice, for other cylindrically symmetric convex bodies, the general method proposed in this 
paper might be even a better choice. 


VI. SPHEROIDS 

Spheroids are cylindrically symmetric ellipsoids (see Fig. Letting a be the semi-axis of the spheroid along the 


(a) (b) 

FIG. 6 . A spheroid is an ellipsoid of revolution. The symmetry axis is here denoted by m; a and b are the ellipsoid’s semi-axes, 
in the direction of m and in the direction orthogonal to m, respectively. This spheroid is said to be prolate (a) if the aspect 
ratio r\ '.— hja is less than unity; it is said to be oblate (b) if 77 > 1 . 

symmetry axis m and b the semi-axis orthogonal to m, we set 

V:=^ (35) 

and call it the aspect ratio of the body. A spheroid with aspect ratio 77 will denoted for short; it is prolate along 
the symmetry axis for 0 < 77 < 1 and oblate for 77 > 1 . Clearly, for 77 = 1 , reduces to a sphere of radius a. 
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Making use of the explicit representation of 5^^ described in Appendix IA3bl we may write the classical Minkowski 
functionals as 


=7ra 2 


Lt 


T]- 


/_! 1 + ( 77 ^ — l)u‘^ 

= 2 iTa^r] J Y^l + (77^ — l)u^du 
V[y^] =yaV =: 14, 


du. 


(36a) 

(36b) 

(36c) 


where Vs has been introduced as a shorthand for the spheroid’s volume. It is often useful to describe how far is 
from a sphere by defining its eccentricity e as 


e := 


Y^l — 77 ^ for 0 ^ 77 ^ 1, 

1 - for 7? ^ 1. 


(37) 


A relevant property of e is that the transformation 771 -^ I/ 77 , which represents the reciprocal inversion of relative 
to its center, changes a prolate spheroid into an oblate spheroid with the same eccentricity. Though neither of the 
functionals (|36]) is invariant under reciprocal inversion of all the ratios 




(38) 


are expected to be so, as such a property should indeed be enjoyed by the ratio of the excluded volume 

of two congruent spheroids, y’l and ^ to their common volumeAs a consequence, all should be functions 

of e alone. The expression for /o was already obtained by Isihara |25| . 


= 2 


l + (l-e2) 


arctanh e 


1 + 


aresm e 
evT^ 


(39) 


which is also known as the Isihara-Ogston-Winzor formula 

The representation for Bn in (|3T]) can appropriately be used to obtain all even-indexed functions /nEE To this end, 
we first record the form taken on a spheroid by the extended Minkowski functional (see Appendix IA8I for more 
details): 


r;^[l+77^ + (l-77y] ^,^ 


Mniy^^] =7ra [ PniO 

^ ^ J -1 + 

M;=27ra Pn (6 

=7ra f ~ ~ y di, 


=27ra‘ 


/>"<«] 


V 


dC 


[ 7?2 + (1 - 7 ?^)^ 2]2 


For 71 = 2, we obtained 

h = 


32? 


3 + (e2+3)(l-e2) 


arctanh ( 


3 - 2e^ + 


4e^ 




: arcsin e 


(40a) 

(40b) 

(40c) 

(40d) 

(40e) 

(40f) 

(41) 


Tjipto-Margo and Evans P credit Holyst and Poniewierski for having proved analytically such an invariance property for uniaxial 
ellipsoids, but we were unable to retrace a convincing analytic proof in |^. Similarly, the extension of this property to biaxial elli^oids 
was established numerically in by a Monte Carlo method. Contrariwise, the explicit analytic formula obtained by Mulder 0, 
for the excluded volume of spheroplatelets allows one to prove that its ratio to the individual spheroplatelet’s volume is invariant under 
reciprocal transformation of the three unequal lengths that characterize these bodies. In any event, as shown in [13, even for spheroids, 
this property does not apply to higher-order virial coefficients. 

Clearly, all odd-indexed fn vanish identically since spheroids are symmetric under central inversion. 
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It is worth noting that this formula coincides with that found by Isihara for oblate spheroids (77 ^ For prolate 
spheroids, Isihara [ 9 | records a result which does not comply with the requirement of /2 being invariant under the 
transformation 77 1-^ I/77. For this reason, we deem it to be incorrect. This should not indeed surprise us, as Isihara’s 
method delivers Bn in the form of two separate power series in e, one for the prolate case and the other for the oblate 
caseEl which need then be resnmmedF^ 

Explicit formulae for both and fe are reproduced in Appendix IA 81 here we shall be contented with showing in 
Fig. ^4, and B2 normalized to Bq as functions of 77 for prolate spheroids (as for oblate spheroids these ratios 

also remain unchanged under the transformation 77 1 /^)- The graphs plotted in Fig. [ 7 | may help deciding how 



FIG. 7. The plots of B 2 (solid line), Ba (dashed line), and Bq (dotted line) normalized to Bo for 0 ^ 77 ^ 1. In the limit 
as 77 ^ 0 (needle-shaped spheroids), the ratios shown here tend to B 2 /B 0 = —5/8 = —0.63, B 4 /B 0 = —9/64 = —0.14, and 
Be/Bo = -65/1024 = -0.06. 

many terms to retain in ([ 4 |) for any given value of 77. 

For completeness, we show in Fig. [8] the graphs of the coefficients Bq and B2, the former of which is normalized to 
8V^, the minimum excluded volume of two congruent spheroids of volume Vs (attained when they are in the parallel 
configuration). 



FIG. 8. (Golor online) (a) The plot of Bo (normalized to 814) for a prolate spheroid. It behaves like 3/3277 as 77 ^ 0. (b) The 
plot of B 2 (normalized to Vs) for a prolate spheroid. It behaves like —157rl32r] as 77 ^ 0. Both plots are easily extended to 
oblate spheroids by preserving their values under the transformation 771 -^ I/ 77 . 


See equations (48)-(50) of Q. 

See equations (29) and (47) of 

A similar discrepancy for /4 is pointed out in Appendix I A 81 below. 
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Over the past few decades, hard ellipsoids have been the object of many studies revisiting the classical Onsager 
theory of hard cylindrical rods. In all these studies, the excluded volume of ellipsoids plays by necessity a central 
role (occasionally, along with some higher-order virial coefficients which have also been computed) 0 More recently, 
a formula was also obtained in [36[ for the excluded volume of two congruent oblate spheroids, elaborating on the 
original method of Isihara That formulaF^ is not directly comparable with ours, as it expresses the excluded 
volume as a power series of trigonometric functions of the angle between the bodies’ symmetry axes, which unlike 
Legendre polynomials is not a system of orthogonal functions. 


VII. CONCLUSIONS 

The major objective of this paper was to express explicitly the excluded volume 14{^ 1 ,^ 2 } of two arbitrary 
cylindrically symmetric, convex bodies and ^2 (with symmetry axes m\ and m 2 ), in terms of shape functionals 
to be evaluated separately for and ^ 2 - We accomplished this task by relating the coefficients that represent 
Ve{^i, ^ 2 } in the basis of Legendre polynomials Pn{mi • m 2 ) to certain anisotropic volume averages which, in 
complete analogy with the classical Minkowski formula for the isotropic average of the excluded volume, were expressed 
in terms of shape functionals that extend Minkowski’s. As demonstrated by the examples of cones and spheroids, 
which we worked out in full details, the extended Minkowski functionals can be evaluated exactly. A large number of 
them might be required to obtain Ve{^i, ^ 2 } at a high degree of accuracy, but the proposed method provides them 
exactly in any desired number. 

As witnessed by the case of cones, one motivation of our study was to explore the role of shape polarity in the 
excluded volume of tapered bodies. It has already been shown that when such congruent bodies and ^2 are convex 
and cylindrically symmetric, Ve{^i^ ^ 2 } attains its minimum in the antiparallel configuration [^. Therefore, one 
could think of assigning a shape dipole d to these bodies by extracting from ^ 2 } the dipolar component, 

Biirii • m 2 , and rewriting it formally as di • Instead, we proved that = 0, thus making elusive the definition 

of any shape dipole for a tapered, cylindrically symmetric, convex body. Clearly, the antipolar property revealed by 
the minimum of Ve{^i, ^ 2 } remains valid, but it can in general be read off from the coefficient ^ 3 , and so properly 
speaking it is an oetupolar effect. 

Cones indeed interested us because they are tapered, but they are not the easiest cylindrically symmetric, convex 
bodies for which one would compute the excluded volume. Perhaps, ellipsoids of revolution might come first in 
anyone’s list. For this reason, we also applied our method to ellipsoids of revolution. Other methods have already 
been devised to compute the excluded volume of these bodies, such as the overlap criteria used in computer simulations 
HIH, or the approximations stipulated in the Gaussian overlap model originally introduced in [45| 0 but with the 
admirable exception of the classical factorized formulae of Isihara [9[ for the first Legendre coefficients Bn and the 
closed form expression for the distance of closest approach for two ellipses in two space dimensions Q 0 no explicit 
analytic representation was known for the excluded volume of ellipsoids of revolution. We hope that we have provided 
one, rooting on geometric grounds the multiplicative structure of Isihara’s formulae and emending some of them. 

Several other applications could be foreseen for our representation formula. In tune again with Onsager’s paper 
[l|, we mention just one: the role of shape in steric interactions of filamentous viruses. This was indeed the original 
motivation of Onsager’s work, which intended to provide a theoretical explanation for the liquid crystalline behavior 
of tobacco mosaic viruses, which were the first to be isolated and purified [l^. An up-to-date review of the recent 
applications of Onsager’s theory to viruses of various elongated shapes can be found in [5T|. We trust that our 
representation formula for the excluded volume could help making the role of viruses’ shape more explicit. 
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appropriate definition of shape dipole for a cylindrically symmetric rigid body, which prompted the study presented 
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20 yyg refer the reader to for a review. Other relevant information can be gathered from the works [33 . [33443| . 

See equation (B12) of |36l| . 

Actually, for selected mi and m 2 on the symmetry axes of the congruent bodies and ^ 2 , one could either orient the vectors di and 
d 2 along mi and m 2 , respectively, or in the opposite directions, provided their orientations are reverted in both bodies. 

An Onsager theory for hard-ellipsoids based on this approximation can be found i n |4al , a paper well aware of the possible inaccuracies 
stemming from the hard-body modification of the simple Gaussian overlap model [4^7 See also [3 for ^ recent review of the Gaussian 
overlap model for hard-ellipsoids. 

Unfortunately, the extension to ellipsoids in three space dimensions of the method that was successful in two dimensions can only be 
performed numerically [i^. 
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Appendix A: Mathematical details 

In this appendix we record for completeness the mathematical details needed to make our development rigorous, 
but which would have hampered our presentation if dispersed in the main body of the paper. We start by recalling 
the essentials of convex body geometry; they are extracted from the wider treatment presented in Appendix A of 
to which the interested reader is referred for a better appreciation of the formalism adopted in this paper. 


1. Essentials of convex body geometry 

A convex body ^ in the three-dimensional space ^ is represented here through the radial mapping u ^ 
which associates to each unit vector u in the unit sphere of S' the point on the boundary of ^ where the 
outward unit normal is precisely Such a representation requires to be mapped univocally onto 9^, which is the 
case whenever ^ belongs to the class of convex bodies with smooth boundaries and strictly positive curvatures. 
Such an assumption is not a true limitation to our development, as is indeed dense in the whole class of 
convex bodies with respect to the Hausdorff metric. Thus, the values attained in by a continuous functional 

defined in can be computed as limits on appropriate approximating sequences of bodies in This property 

is for example exploited in Se c. lA 4 ^ below to compute the contribution of a sharp ridge to the extended Minkowski 
functionals introduced in Sec. lIVFi 

Figure [9] illustrates our representation of SS through its radial mapping r. It also shows that the unit outward 



FIG. 9. (Color online) Sketch that describes how the radial mapping r assigns to a unit vector v of the translation that 
brings o into the point on where v is the unit outward normal to d^. The existence of such a mapping is guaranteed by 
the assumption that ^ belongs to . 

normal n to 9^, which by construction at the point r{iy) coincides with i/, can also be regarded as a field on 
Its surface gradient Vgn is the curvature tensor and can be represented as 

Vgn = cTiei O ei + 0-262 O 62, (AI) 

where the positive scalars ai and 0-2 are the principal curvatures of 9^, and the orthogonal unit vectors ei and 62 , 
both tangent to 9^, designates the principal directions of curvature. In this paper, fully devoted to cylindrically 
symmetric bodies, we have conventionally taken ei along the local meridian, so that ei, n, and the symmetry axis m 
of ^ are everywhere in one and the same plane (possibly varying with the point selected on d^). Figure [TOl shows 
the geometric situation envisaged here. 



FIG. 10. (Color online) Cross section of a cylindrically symmetric body ^ through a plane containing its axis of symmetry m. 
Both the outer unit normal n and the principal direction of curvature ei along the local meridian are on this plane. 
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The very same property makes it possible to arrive at the expressions for the extended M and S functionals of a cone listed in 
Sec.E 
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The mean curvature H and the Gaussian curvature K are defined in terms of the principal curvatures as 

H \= \ (cTi + CF2) and K := (Ji(J2. (A2) 

The former can also be expressed as 

= tr Vgn = ^ divs n, (A3) 

where tr is the trace operator and divg denotes the surface divergence. Similarly, letting A* denote the adjugate of a 
second-rank tensor we also have that 

(Vgn)* = Kn (8) n = Ku 0 v. (A4) 

The surface gradient VsT* of the radial mapping r over has an expression similar to dsn), 


where 


Vgr = piei (8) ei + ^262 (8) 62, 


Pi := — and p 2 := — 

(Ji (72 


are the principal radii of curvature of d^. In complete analogy with dSij, we have that 

(Vgr)* = 

whence it follows that the surface dilation ratio induced by the mapping r that sends onto is given 

^ = Kv.ov, = 1 

Putting together (IA5II . (IA6I) . (I A 31) . and (IA2I) . we can also write 

divs r = pi + p 2 = ^ divs n. 

K 


(AS) 

(A6) 

(A7) 

I 

(A8) 

(A9) 


In the following, we shall also denote by x the position vector on d^. Formally, the fields u and n are related through 
X by the relations 


u = n{x) and x = r{u). 

A theorem that we have often used in this paper is the surface-divergence theoremF^ It says that 


/ divs uda{n) = / (divg n)u • nda{n) = 2 / Hu • nda{n), 

Jy’ Jy’ Jy’ 


(AIO) 

(All) 


for any continuously differentiable field u defined on a closed smooth surface ^ with unit outer normal n and mean 
curvature H. 

Three continuous functionals defined on the whole class of convex bodies were introduced by Minkowski. They 
are the total mean curvature M, the surface area S, and the volume V. For a body ^ G they are defined and 

represented as follows 0 


M[^] := [ Hda{n) = 

Jd3S 

S[M\ := f da{n) = f —da{u) 
Jd3S K 


/ r • uda{u)^ 


[ V • (Vs'r)*i/da(i/), 


(Al2a) 

(AI2b) 


A* is characterized by requiring that A*(u X v) = Au X An, for all vectors u and v, where X denotes the cross product of vectors (see 
also Sec. 2.11 of B). 

See also Sec. 5.2 ofl^. 

See also Sec. 5.2.3 o f[^ . 

See Appendix A of |1^ . for more details 
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V[^] ^ f n'Xda{n) = ^ f r • iy^da{u) = ^ f {u • r)u ■ {V^rYuda{u). (A12c) 

3 Jam 3 jg2 A 3 Jg2 

As shown in greater details in HI, one of the advantages of representing a body ^ in through its radial 
mapping r is that the Minkowski sum +^2 of two bodies, and ^ 2 , represented by the radial mappings ri and 
r 2 , respectively, is represented by the radial mapping ri 2 = ri + r 2 . Correspondingly, the fundamental functionals in 
(|A12p evaluated on the Minkoski sum of two bodies, and ^ 2 ? of are delivered 

M[S§i + S§ 2 \ = M[^i] + M[^ 2 ], (A13a) 


S[3§i + ^ 2 ] = S[^^i] + S[^2] + 


f [sin2 

L 


(AVf ’ + P2^P2^) + 




da{v), (A13b) 


V[S§i + S§2] = V[S§i] + V[S§2] 






ri 


K(^) 


1 / ■r2 


kw 


da{u) 


+ 3 ■ '^1 + ^ ■ '^2) [sin^ (f> [pY’pY’ + P2 ^p‘i^) + cos^ (f> [pY'p'2' + P2 ’pY’)\ da{v), 


(1) (2) , (1) (2)\ 


(A14) 


where and are the principal radii of curvature of 9^i, and p^ 2 ^ are those of 9^2, 4> ^ [0 ? ^tt] is the 
angle of the rotation about v that brings the pair of principal curvature directions ^ 2 ^^) of into the pair of 

principal curvature directions of body ^ 2 , and = {p^i^ p^ 2 ^)~^ ^ = {p^i^ are the Gaussian 

curvatures of d3Si and 9=^2, respectively. 

We finally remark that for a body ^ G represented by the radial mapping r{u)^ the central inverse (relative 
to the same origin o) is represented by the radial mapping r* defined by 


r*(i>') := —r{—iy). (A15) 

As a result, if ri and r 2 are the radial mappings representing the bodies and ^2 in the body 

whose volume, by Mulder’s identity ©, is the excluded volume Ve{^i^ ^ 2 } of the pair (^ 1 ,^ 2 ), is represented by 

the radial mapping 


re(i^) := ri{u) - r 2 {-v'). (A16) 

It is not difficult to show with aid of (jA15p that the shape functionals defined in (|24|) for ^ G are invariant under 
the transformation ^ i-G- . 

Another conse quen ce of Mulder’s identity is that (jA14[) bears a close resemblance to Wertheim’s representation for 
Mayer’s function An important difference, however, between our method and Wertheim’s is that the expansion 

in ^ with coefficients Bn as in (|3T|) does not result from an expansion of the integrand in the last integral of (jA14p , 
thus avoiding the ambiguities acknowledged in [53| HI 


2. Anisotropic volume averages 


The anisotropic volume averages (PnV) [=^ 1 ,^ 2 ] are defined in ([T3j). The first average {PiV) [^i,=^ 2 ] has been 
computed in Sec. uni here we compute all others. The method employed will be the same as in Sec. uni but to make 
it effective we need to replace (dD with the more general addition formula (see Sec. 18.18.9 of HU), 






30 See (A25), (A43), and (A4^of O. 

3^ See, in particular, (36) of |^, which in an incomplete form was also referred to as the convolution decomposition of Mayer’s function 
by Rosenfeld and later re-established in [30| in its complete form, equivalent to Wertheim’s original equation. 

32 Compare, for example, (64) and (68) of |^. 
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where is the Jacobi polynomial of degree n and indices Jacobi polynomials are defined in the interval 

[—1,1] and are orthogonal relative to the weight function w{x) = (1 — x)^{l-\-xY. They enjoy the symmetry property 
p^a,/3)(_^) _ and can be represented as finite sums (see Sec. 18.5.8 of [ 2 l|), 




(A18) 


The first three Jacobi polynomials that interest us are 


= 1, = 3a:;, P 2 ‘^’‘^’{x) = 7x^ — 1 


?( 2 , 2 ) 


( 2 , 2 )/ 


(A19) 


With the aid of (1 1811 and (|A17I) . we establish the identity, 


T P„(mi • m 2 ) sin^ (p ^ + cos^ p {p^^ p2^ + # 

= P„(mi ■ u)Pn{m2 ■ v) 

“ “ ("^1 ■ - («^2 • iyf]Pi^7.l\mi ■ iy)Pi^4\^2 ■ V) (^Pi^ - P2^^ (pi^ - P2^) , 


(A 20 ) 


where, as stipulated above, the principal directions of curvatures and for bodies and ^ 2 , respectively, to 

which the principal radii of curvature and are correspondingly associated, lie orderly on the planes (mi,n) 
and (m 2 ,n). Use of (pT|) . ()A17j) . and (|A 2 Q[) in (jA14[) leads us to 


(PnV) 


h,^2] = - 


\ (/ Pn{m2 ■ v) 


i^(2) 


/ {u ■ ri)Pn{mi ■ u)da{u) + {{u ■ r 2 )Pn{m 2 ■ v))^ / 

V JS2 JS2 


Pnjmi ■ v) 

kw 


+ 


i (Pn{m2 ■ v) {pf^ + P 2 ^))^ jjp ■ ri)Pn{mi ■ V (pP + da{u) 

+ ^ ((*^ • r2)Pn{m2 ■ V Pnirni ■ v) ^ da{v) 

1 (n - 2 )!(n + 2 )! / , 2173 ( 2 , 2 )/ ,/ ( 2 ) ( 2 )\\ 

' [1 - (m 2 • V) ]P,)^2 («^2 • V) -P2 ) y 


6 (4n!)^ 

X - (mi ■ ■ ri)P^^_lf{Tni ■ u) (pP - da(i/) 

“ ( 4 ^ 02 ^^^' ■ ^^^^n-l\'m2 ■ I') (pf^ - P2^))^ 

X - (mi ■ uf]Pp^{mi ■ u) (pb - da{u). 


(A 21 ) 


To accomplish our task we need now compute all the integrals featuring in (jA2ip . To make this easier, it is expedient 
to realize that they result from parameterizing some general shape functionals through the radial mappings ri and 
r 2 of bodies and ^ 2 - For a cylindrically symmetric body ^ G by (jA 8 p . we see that 


J§2 A Jq 




Pn{Tn • n)da{n). 


(A 22 a) 
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Similarly, also by use of (IMl), (IMl, and (IMl, we easily arrive at 

/ {u - r)Pn{rn - u)da(y) = / {u ■ x)Pn{m ' n)Kda{n)^ 

/ Pn{m ■ u){pi P p 2 )da{i/) = 2 / Pn{m • n)Hda{n), 

Jd^ 

/ {ly ' r)Pn{m ‘ v'){pi P p2)da{v') = 2 / {n ' x)Pn{m ' n)Hda{n) 
^S2 Jd^ 

( [1 - (m ■ ■ r)Pp\m ■ i/)(pi - p 2 )da{v) =-2 [ [1 - (m ■ n)2](n ■ x)P, 


( 2 , 2 ) 

n—2 


(A22b) 

(A22c) 

(A22d) 


(m • n)\{cFi - G2)da{n), 


(A22e) 

f [1 - (m • i/)^]P^^ 2 ^(m • - p2)da(y) = -2 / [1 - (m • n)^]P^^ 2 ^(m • n)\{(Ji - cF2)da{n) (A22f) 

In formulae (|A22p we readily recognize the shape functionals defined in (|24|) . With the aid of these definitions, we 
give (jA2ip the form (|25|) used in the main text. 


3. Generating curve 

Here we represent the boundary of a cylindrically symmetric convex body ^ as generated by the 27r-rotation 
of a plane curve, y{s) = r{s)ex — a{s)ez, parameterized in the generic scalar s (see Fig. [TT]) . Identifying m with the 



FIG. 11. (Color online) The plane curve y, which generates by a 27r-rotation about m, is parameterized in the generic 
scalar s as y{s) = a3(s,0), where x{s,(f) is given by ()A23D . The origin o is taken to coincide with the uppermost pole. 

unit vector Cz of a Cartesian frame (e^^, Cy^Cz)^ we can then represent as the surface 

x{s, ip) = r{s)er — a{s)ez, (A23) 


where 


Cr = COS ipCx P sin ip By 


(A24a) 


is the radial unit vector and 


= —fillip Bx P ip By (A24b) 

is the associate orthogonal unit vector in the plane (e^^. By). 

By letting s and ip depend on a parameter t, we obtain a trajectory t ^{t) := x{s{t),ip{t)) on d^. It follows 
from (|A23p that 


i = sy/ + ipvB^, (A25) 

where a prime ' denotes differentiation with respect to s, a superimposed dot denotes differentiation with respect to 
t, and 


pBr — a' Bz 

\/P^ + 


(A26) 
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is the unit tangent vector to cc(', (/:?), for given ip. From (|A26p and (|A25p . we easily arrive at both the unit outward 
normal to 


n 


a' Br + t'Bz 


and the surface area element 


da{n) = r\/ + a'‘^dsdp. 

By further differentiating n along the trajectory ^(t), we obtain that 


s{a"r' — a'r") pa' 


n = 


Since h = (Vgri)^ and, by (jA25p . 


s = 


+ a'^ 






and p = 


for ^ is arbitrary, we conclude that 


(A27) 


(A28) 


(A29) 


(A30) 


V.n = 


rt ( 8 ) t + ■ 


(r'2 + " ' r\/r'^ + 

whence we read off at once the principal curvatures of d^. 


(A31) 


a. Cones 

Figure [12] depicts the generating curve for a circular cone with vertex in the origin o, semi-amplitude a, radius 
R and height h, which are related to the slant height L through the equations 



FIG. 12. (Color online) The generating curve of a circular cone with semi-amplitude o, radius R and height /i, the two 
latter related to the slant height L as in (IA32D . The parameter s here represents the arc-length along the slant side of the cone. 

R = L sin (a, h = L cos a. (A32) 

The functions r{s) and a{s) featuring in (jA23p are correspondingly given by 

r{s) = s sin a and a(s) = scos(a, (A33) 

where now s has been chosen as the arc-length along the slant height of the cone; it follows from (jA3ip that 

, cot a / * . X 

(Ji = 0 and <72 =-. (A34) 
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FIG. 13. (Color online) The half-ellipse with semi-axes a and b whose 27r-rotation about the symmetry axis m — Cx generates 
a spheroid with aspect ratio 77 = h/a. The parameter 0 featuring in (|A35D designates the angle between y{0) and Cz- The 
origin o has been chosen in the center of the spheroid. 


b. Spheroids 


The generating curve for a spheroid is illustrated in Fig. [131 is a half-ellipse with semi-axes a and 6 , along 
Cz and Cx respectively, and centered in the origin o. Letting the parameter s be the angle 0 ranging in [0, tt] and 
depicted in Fig. [131 the functions a{s) and r{s) in (|A23p are now written a® 


a{0) = —acosO and r{0) = 6 sin^. 
By use of (jA27p and (|A3ip . we readily arrive at 
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a [1 

+ 

(^2 _ 

1) cos^ 


1 
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av \J\ + (jf' - 1) cos^ Q' 


(A35) 

(A36a) 

(A36b) 

(A36c) 

(A36d) 


where rj := b/a is the spheroid’s aspect ratio. 

4. Extended M and S functionals of a circular ridge 

Here we apply the formalism presented in Sec. lA3l to compute the extended M and S functionals defined in Sec. m 
for a circular ridge 91 of radius i?, where neither H nor K are defined. To this end, we replace 91 with a toroidal 
approximation 91^ with equatorial radius R and meridian radius 5 , whose outer unit normal n spans the sector in 
which the angle 0 that it makes with the symmetry axis m ranges in the interval [ 6 >i, 6 ^ 2 ]. To afford a greater generality 
(and in view of our application to cones in Sec. [V] above), we choose the origin o on the symmetry axis at the generic 
distance h from the ridge’s plane (see Fig. n. Our strategy will be to compute the extended M and S functionals 
on 91^ and then take the limit as 5 ^ 0+. The functions r{s) and a{s) introduced in Sec. lA 31 which here describe 91^ 
are 


r{s) = R£ sin 0{s), a{s) = h— £ cos 0{s), (A37) 

where 0 and s are related through s — si = £{0{s) — 6 >i), with si an arbitrary constant. It easily follows from (jA3ip 
and (|A37p that 


1 

5 


sin^ 


(^2 = 


R^ esinO 


(A38) 


To avoid typographical clutter, we are guilty of using the same symbol for both the function a{6) and the scaling semi-axis of the 
generating half-ellipse. 
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FIG. 14. (Color online) (a) Circular ridge D1 of radius R in the plane orthogonal to the symmetry axis m at the distance h 
from the origin o. The unit outward normal n makes the angle 0i with m on one side and angle O 2 on the other side, (b) The 
sharp corner of 91 is rounded off in a toroidal surface with meridian radius £. 


Moreover, (jA28[) yields 


da{n) = {R e sin 0)ed0 dip. 


(A39) 


Using (jA38[) and (jA39|) in the definitions of the extended M and S functionals in (El, and then taking the limit 
as 5 0+, we arrive at the following expressions: 


r02 

Mn[^]=7rR Pn (cos (9)d<9, (A40a) 

Joi 

r02 

M^[91] = 27r / {RsinO — h cos0) sinOPn{cos 0)d0, (A40b) 

Joi 

p02 

M"[9i]=7ri?/ sm^dP^^^l\cosd)dd, (A40c) 

Joi 

Snm=0, (A40d) 

r02 

5'^[91]=7rP / {RsinO — h cos0)Pn{cos 0)d0, (A40e) 

Joi 

r02 

[^] = 0 — h COS 0) sin^ ^Pn -^2 0)d0. (A40f) 

Joi 


a. Extended M and S funetionals for a disk 


Formulae ()A4Qp are instrumental to obtaining the explicit expressions for the extended M and S functionals of 
a disk ^ of radius R. As before, we start by replacing ^ with an approximating rounded body, the spherodisk 
defined as the Minkowski sum of ^ and a ball of radius 5 and center coincident with the center of Figure [15] 
illustrates both and the generating curve of its boundary. The extended M and S functionals for & will be 
obtained by taking the limit as 5 ^ 0+ in those computed for consists of two fiat parallel disks, for which 

both principal curvatures vanish, and the toroidal approximation 91e of the circular rim 91 of for which the angles 
0i and O 2 in Fig.[T4]are = 0 and O 2 = tt, respectively. Apart from the limit as 5 ^ 0^ of Snl^s]^ which is immediate 
to compute, for all other functionals this limit follows directly from (|A4Qp by setting h = 0 and choosing 0i and O 2 as 
















22 


m 

( e) 


FIG. 15. (Color online) (a) Spherodisk Sig defined as the Minkowski sum of the disk & and a ball of radius e and same 
center o as (b) The generating curve of The symmetry axis m is orthogonal to 


above. We thus arrive at 


Mn[&] =7rR [ Pn(cos(9)d<9, 

Jo 

(A41a) 

= 27rP f Pn{cosO) s'm^ OdO, 

Jo 

(A41b) 

M;' [^] =7rR [ P^^df (cos 0) sin^ OdO, 

Jo 

(A41c) 

=7ri?2 (P„(l)+P„(-1)), 

(A41d) 

S'j^[S’] = ttR^ I Pn{cosd) sin9d9, 

Jo 

(A41e) 

5"[^] = 7rP2 f p^^2) 

Jo 

(A41f) 


In particular, it follows from (jA4ip that all extended M and S functionals with an odd index n vanish for a disk. 



5. Invariance under translations 


The anisotropic volume averages for which we found in (|25]) an explicit representation in terms of the extended 
Minkowski functionals are clearly invariant under the full Euclidean group comprising both translations and rotations. 
On the other hand, as clearly shown by equations (|24|) . while all extended M and S functionals are invariant under 
rotations, those that also appear to be invariant under translations are only M^, M^', and 5'^- 

are expressed as integrals over the boundary of the body ^ of fields that depend explicitly on the origin 
o through the position vector x. Here we shall show that, despite all appearances, is indeed invariant under 
translations, whereas both S'^ and 5'^' are not. For the latter two, we shall also give explicit formulae that describe 
how they are affected by a translation. Of course, the combination of these functionals in (|2^ must be translation- 
invariant. We shall exploit this fact in Sec. lA 61 below to show that functionals M'^ and are not independent, a 
conclusion which would be hard to reach by direct comparison of their definitions. 

Translating a body ^ by the vector a is formally equivalent to taking the Minkowski sum ^ + a of ^ and the point 
in space identified by a. Moreover, since all extended M and S functionals are invariant under rotations, computed 
on ^ + a for any given they are isotropic functions of a. It readily follows from (j24b[) that 

M;[^ + a] = M;[^] -b a • / Pn{m- n)Knda{n). (A42) 

Jam 

The integral on the right side of (|A42p is an isotropic vector-valued function of m; as such, by the Cauchy theorem 
on isotropic vector-valued functions, it must be proportional to m. Thus, (|A42p becomes 

+ a] = M^[^]-f-a • m / {rn • ri)Pn{rri • ri)Kda{ri). (A43) 

Jd3S 

For ^ G by use of (jASp . we see that 


f {rn ■ n)Pn{rn ' ri)Kda{n) = f Pi{m ■ iy)Pn{m ■ iy)da{iy) = 27r f Pi{x)Pn{x)dx = 0 V n ^ 2, (A44) 

Jd^ J-1 
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where the last equality follows from the orthogonality of Legendre polynomials. Since we have already proved in 
Sec, mil that Mi[^] vanishes identically for all ^ G by (jA44p we conclude that all functionals are invariant 
under translations. 

This is not the case for both S'^ and S'J^. Reasoning precisely as above and making use of the recurrence relation^ 


n + 1 


xPnix) - + + 2n + 


xp: 


( 2 , 2 ) 


n—2 


^ (n-l)(n + 3) p( 2 , 2 ). s n ^( 2 , 2 ). . 
^ ’ (n + l)(2 n + 1 )^ ^ ^ 2n + 1 ^ 

the latter valid for n ^ 2 and with the postulation that = 0, we arrive at 

(n-l)(n + 3) 


51 


■ al = 5: 


S'l\^% + a] = S'l 


+ a - m 


^ (ti + 1 ) (2n + 1) 

the latter valid for n ^ 2 and with the postulation that = 0 


M"+i[^] + 


277/ T 1 


K-i[ 


(A45a) 

(A45b) 


(A46a) 

(A46b) 


6. Reduction formulae 

Here we take advantage of the general formulae (jA46p just established and of the specific expressions for the 
extended M and S functionals obtained in (jA4ip to show that each functional reduces to Mn and to substantiate 
our conjecture that so should equally do each M^. 


a. M'n reduced to Mn 

By requiring that the anisotropic volume averages, as expressed by (ESI, be invariant under translations for all 
bodies and ^ 2 , a laborious but easy computation relying on (jA46p and the translation-invariance of shows 
that 


where the coefficients must obey the recurrence equation 

{n + 3){n + 2) 


16(77/ + 1)77/ 


whose explicit solution is 


^n+l^n ~ 


477/ 


77/ -|- 2 


(A47) 


(A48) 


(A49) 


Combining (jA49p with (jA47p . we arrive immediately at (j27j). 


b. M' reduced to Mn 


Inspired by (jA49p . we computed the ratio a'^ of to Mn[S^] for a disk interpolating with the aid of (jA41aP 

and (jA41bp the values of obtained for a number of indices tT/, we concluded that 


(tT/ — l){n + 1) 


V 77/ > 2, 


(A50) 


See, for example, Sec. 18.9.1 of [2ll| . 
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whence (|28]) follows at once. Although we could not establish (|28|) on a firmer basis, we checked by use of (j33ap and 
(|33bp and of (|4Qap and (|4Qbp that it is valid for a large number of indices n when ^ is taken to be either a cone 
or a spheroid for all values of the semi-amplitude a and of the aspect ratio r]. We are aware that and pSP 
have a completely different standing, as the former has been proved rigorously, whereas the latter is only conjectured. 
Most of our development in the main body of the paper relies neither on nor on pS]) . What does depend on pSP 
is only the possibility of giving compact factorized formulas for the coefficients Bn as those listed in (jA5ip and (jA54p 
for cones and spheroids, respectively, both of which are expected to obey pSp . 


7. Legendre coefficients for the excluded volume of cones 

Letting and ^2 be two congruent circular cones, and ^ 2 ^ with semi-amplitude a, with the aid of p^ . pip . 
and p3P we arrived at the following explicit formulae for the first eight Legendre coefficients Bn plotted in Figs. [3] 


andlU^b) as functions of a: 

Bo = sin^acosa + sin a + a) sin a + cos a] (1 + sin a), (A51a) 

Bx = 0, (A51b) 

B 2 = sin(a(2asina + Trsina -1- 2 cos a — 6 cos^a)(3 cos^a — 2 — 2 sin a), (A51c) 

Bo = sin^acos^a, (A51d) 

^4 = ~ 2 ^^^^ sin (a(37r sin a 6a sin a 6 cos a — 130 cos^a + 140 cos^a) 

X (35cos^(a — 40cos^(a + 8 + Ssina), (A51e) 

Bo = sin^acos^a(27cos^(a — 20)(9cos^a — 8), (A51f) 

Bq = 555^36 sin a{b7r sin a + 10a sin a + 10 cos a — 686 cos^a + 1876 cos^a — 1232 cos^a) 

X (231 cos^a — 378 cos^a + 168 cos^a — 16 — 16 sin a), (A51g) 

Bj = sm^acos^a{2S0 — 924cos^(a + 715 cos^ a)(143 cos^a — 198cos^a + 72), (A51h) 


where L is the cone’s slant height. They are recorded here both for completeness and as an illustration of the method 
proposed in this paper. 


8. Legendre coefficients for the excluded volume of spheroids 


To obtain the coefficients Bn that express the excluded volume of congruent spheroids as a series of Legendre 
polynomials, we computed the extended Minkowski functionals in pTp for the generating curve described by (jA35p . 
Use of (|A36p in (|24aD gave 


Mn[y'^]=TTa 


/■ 


TfU 




-1- 1 1 du^ 


where we have set u := cos^. The change of variables 

i •= 


rju 


y/l + {rf - l)u^ 

then led us from (jA52p to (j4Qap . The other formulae in pOp were obtained in precisely the same way. 
Along with the expression for /2 recorded in (|4Tp . we also obtained 
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In Fig. Uni using (|A54aP and (|3^ we plotted the ratio = / 4//0 against 77 . It is there also contrasted against 



FIG. 16. (Color online) The ratio B^jB^ — / 4//0 is plotted against 77 both according to the expressions in (IA54a|) and (|39l) 
(solid line) and to formula (60) of (dashed line). 

the function obtained for this ratio by Isihara [9| (only for the prolate case). The two graphs fail to coincide, even 
dramatically so, away from 77 = 1. In particular, we estimate that 

lim —^ = 9, (A55) 

77^0 d(I) ^ ^ 

where b!^ is B^ as delivered by equation (60) of [9[. 


Appendix B: Shape-reconstruction method 

In this appendix we describe the method adopted for reconstructing the boundary of the excluded body 
for two congruent circular cones and semi-amplitude a and slant height L, hereafter simply denoted for 

short. More precisely, the method reconstructs a triangular surface mesh that, depending on a fundamental parameter 
to be described, approximates at any degree of precision. From the surface mesh, the approximate value of the 
excluded volume can be computed immediately. 

The method adopted for this task is a pipeline of two algorithmic components: 

1 . an online vector quantization algorithm that includes a generator of random point samples from and which 
produces a configuration of reference vectors W ; 

2 . a surface reconstruction algorithm that produces from W the triangulated surface mesh that represents an 
approximation to 

The method described is similar to that in 0 - In particular, the random generator of point samples from d^e is 
essentialy the same. In that context, however, all target surfaces d^e were generated from sphero-cones and could 
be assumed to be smooth, so that the reconstruction process could be embedded into step 1 above via the SOAM 
algorithm By contrast, in the case of cones considered here, the presence of ridges and cusps in forces 

adopting a different strategy. In the rest of this appendix, the main aspects of this new strategy are discussed in 
detail. 

Upon comparing the method described here with others computing the densest packing of particles of arbitrary 
shape [iH, we heed in passing that our method determines directly the surface bounding the excluded region, with 
arbitrary degree of precision and in one run per pose, whereas those other methods typically require repeated Monte 
Carlo simulations 1^. 


1. Sampling the surface boundary 

Random point samples from d^e can be generated with a procedure based on equation (jA16p . reproduced here for 
convenience: 


re(iz) = ri(iz) - r 2 (- 2 ^). 
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Here reaches a point on d^e and ri and r 2 , in this specific case, designate points on and respectively. 

Random points on can be obtained either by a generating a random vector reaching a point on and then 
finding a vector to a point on d ^2 opposite normal — i/ or by reverting this very procedure: the sum of the 

vectors thus obtained will belong to d^e. 

The main difficulty in implementing such a random generator is to guarantee positive sampling probability almost 
everywhere on that is, apart from subsets of zero area measure. On all smooth components of a circular cone, in 
fact, the Gaussian curvature K vanishes and this means that in general a normal vector ly does not identify uniquely 
one point on the cone’s surface. Furthermore, the Minkowski sum of two straight lines on the boundary of each cone 
can result in a surface patch with positive area measure on despite the fact that each line has zero area measure 
and thus no chances of being sampled, unless specific provisions are introduced. Appendix B in [T^ describes how 
these problems can be circumvented in actual computations. 

Although the requirement of positive sampling probability almost everywhere can be enforced in practice, no known 
method guarantees uniform sampling probability over As shown in Fig. Era), the overall sampling obtained 

with the chosen random point generation method is indeed non-uniform. 



FIG. 17. (Color online) (a) A set of 30 K random points generated from d^e with the method adopted here: the sampling 
of the surface is clearly non-uniform; (b) the final configuration of lOK reference vectors W produced by the adaptive NG 
algorithm is more uniformly distributed; (c) and (d) from W, the ball-pivoting algorithm reconstructs the surface boundary 
d^e with no human intervention. 


2. Vector quantization: adaptive neural gas 

Many well-known algorithms for surface reconstruction work considerably better when the input point cloud is as 
close as possible to a uniform sample of the target surface and are often hampered when this is not the case@ Apart 
from greater time complexity, these difficulties can lead in practice to the need for accurate verification of results and 
possibly to manual post-processing, to correct imperfections. 

The intended purpose of a vector quantization algorithm in this context is to obtain both an improvement in the 
uniformity of sampling and a quantitative reduction in the number of points to be used for surface reconstruction. 
The algorithm of choice is an adaptive variant of the neural gas (NG) algorithm [59| and works as follows: 

1. initialize W with a pre-defined number k of reference vectors Wi positioned at random on 

2. generate a random point p from 

3. find the nearest reference vector in FF, i.e. Wi := argmin^^^.^vv ||p “ '^ilh 

4. if ||p — Wi\\ < r, where r is a fixed threshold, adapt all reference vectors in W by 

Awi = e • hx{ki{p)){p - Wi) 

where ki{p) := ^{wj : ||p — mj|| < ||p — mi||} denotes cardinality), 5 > 0 is a real parameter and 

ho{n) := (ion and hx{n) := for A > 0; 


Known methods for uniform sampling presuppose knowledge of the surface’s analytic description plus further specific conditions [57|. 
More precisely, the relevant requisite in this respect is that the point sampling should be at least locally uniform |58l| . 
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5. otherwise, if p is farther away from add a new reference vector pioW] 

6 . unless a maximum number of iterations T has been reached, return to step 2. 

As evident from step 5, this algorithm is adaptive in the number of reference vectors in W; in particular, this 
means that the level of refinement of the sampling of d^e provided by W can be controlled through the value of the 
fundamental threshold r. 

In [m it is proven that, when the value of the constant 5 tends to 0 as the iterations progress, the NG algorithm 
performs a stochastic gradient descent towards a (local) minimum of an overall cost function and that its configuration 
tends to obey the power law 


p{w) (X P{wy with 7 -=^^: 

where d is the dimension of the input space being sampled, that is d = 2 in this case. Here p{w) is the density of 
reference vectors inWatw and P{w) is the sampling probability. Since the exponent 7 is smaller than 1 , the overall 
configuration of W tends to be closer to uniformity than the sampling probability P. This effect is clearly visible in 

Fig.HZfb). 


3. Surface reconstruction 

With proper parameter settings (see below), the reconstruction of a triangular mesh from the final configuration 
W produced by the adaptive NG algorithm poses no particular problem and could be performed in full automation. 
In this work we used the ball-pivoting algorithm [b^ which joins in a triangular face any three vectors in W whose 
ends are touched by a ball of a given radius r that does not contain any other vector’s end from the same set. One 
example of the results of this procedure is shown in Figs. HHc) and (d). Further examples are shown in Fig. [181 which 
contains a gallery of shapes produced with the method described above. 


/ ^ / 

(a) (b) (c) 


m c> u 

(d) (e) (f) 



(g) (h) (i) 

FIG. 18. Gallery of reconstructed boundaries d^e for pairs of congruent circular cones of semi-amplitude a. Rows correspond 
to values of a equal to and I|7r, respectively, while columns correspond to values of the angle'd between the symmetry 

axes mi and m 2 equal to 0, f and tt, respectively. All figures are in the same scale and frame of reference. 
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4. Implementation and benchmark 

The adaptive NG vector quantization algorithm, together with the generator of random points from 83^^^ has been 
implemented in Java. In order to speed the execution up, the algorithm has been converted to a multi-threaded 
version suitable for multi-core computers, along the lines described in 0. For surface reconstruction, we used the 
implementation of the ball-pivoting algorithm included in the Meshlab open-source tool [HU . 

The overall method for shape reconstruction was validated using Minkowski’s formula for isotropic volume average 
m together with the cone-specific functionals (|3^ . For benchmarking, pairs of congruent circular cones and ^2 
having slant height L and semi-amplitude a varying from ^ to ^tt with step ^ were considered. For each such pair, 
the value of V[^e] was computed for angles'd between the two symmetry axes mi and m 2 varying from 0 to tt with 
step the isotropic average of the resulting sequence of volumes was then computed and compared with the exact 
value of (y)[^e]- The fundamental threshold r, which governs the density of reference vectors in W with respect 
to 9^6 5 was determined empirically with the objective of having a difference lesser than 0.02% between the exact 
value of each isotropic average and the corresponding value computed numerically. A value r = -^L was found to be 
adequate (see also the comparative plots in Figs. lU^a) and (b)). Also the value of T = 120 M maximum equivalent 
iterations of the NG algorithm was determined empirically. In the actual experiments, the execution was split into 4 
concurrent threads, each processing in multi-signal mode (see 0 ) 250 random points per iteration. Being dependent 
on the area of the number of reference vectors in the final configurations of W varied greatly, from 3,592 to 

41,689. 

All numerical experiments were run on a workstation based on an Intel® Xeon® GPU E3-1240 v3, 3.4 GHz GPU 
with 8 GB of RAM. As for computing times, the most demanding part of the method is running the T = 120 M 
equivalent iterations of the adaptive NG algorithm. For each pair of cones and for each pose, with the precision 
required, this computation took on average about 4,254 seconds (i.e. about 71 minutes) to complete. 
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